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In this paper, we demonstrate the use of physical models to evaluate the production of ^^Ar and 
^"Ar underground. Considering both cosmogenic ^^Ar production and radiogenic ""'Ar production 
in situ and from external sources, we can derive the ratio of "^^Ar to ''"Ar in underground sources. 
We show for the first time that the ^® Ar production underground is dominated by stopping negative 
muon capture on '^^K and (a, n) induced subsequent ''^K(n,p)^®Ar reactions. The production of 
''^Ar is shown as a function of depth. We demonstrate that argon depleted in ^^Ar can be obtained 
only if the depth of the underground resources is greater than 500 m.w.e. below the surface. 
Stopping negative muon capture on '^^K dominates over radiogenic production at depths of less 
than 2000 m.w.e., and that production by muon-induced neutrons is subdominant at any depth. 
The depletion factor depends strongly on both radioactivity level and potassium content in the rock. 
We measure the radioactivity concentration and potassium concentration in the rock for a potential 
site of an underground argon source in South Dakota. Depending on the probability of ^^Ar and 
''"Ar produced underground being dissolved in the water, the upper limit of the concentration of 
■'^Ar in the underground water at this site is estimated to be in a range of a factor of 1.6 to 155 less 
than the "^''Ar concentration in the atmosphere. The calculation tools presented in this paper are 
also critical to the dating method with '^^Ar. 

PACS numbers: 95.35.+d, ll.lO.Lm, 29.40.Mc 



I. INTRODUCTION 

Argon is expected to be an excellent target for direct 
dark matter detection experiments searching for Weakly 
Interacting Massive Particles (WIMPs) (HQ- However, 
natural argon collected from the atmosphere contains 
■^^Ar and ^^Ar that are radioactive isotopes. The ratio of 
^''Ar to ''^Ar is measured to be 8.1 x 10^^'^g/g [Sj-^] for at- 
mospheric argon. Such a concentration results in 1 Bq/kg 
of specific activity of ^^Ar (Q = 565 keV, ti/2= 269 y|ei). 
Existing experimental limits on the ^^Ar concentration in 
the atmosphere are less than 10^^^ '^^Ar atoms per '"^Ar 
atom [3-0|: which makes about 0.018 Hz/kg of specific 
activity of ''^Ar (Q = 3520 keV and ti/2 = 33 y). The 
production of '^^Ar in the atmosphere is dominated by 
the ^°Ar(n, 2n) process [s*] with a reaction threshold of 
1011.86 keV. ^^Ar is produced in the atmosphere through 
""^Ar -t-n ^^Ar and ^^Ar+n — )-*^Ar. The neutrons in 
the above reactions are from cosmic rays. 

Recently developed argon-based dark matter detectors 
measure scintillation light induced by low-energy nuclear 
recoils due to elastic scattering of WIMPs. Argon has two 
distinct mechanisms for the emission of scintillation light: 
prompt light (6 ns) and delayed light (1.5 /is) 10]. The 
ratio of the prompt light to total light has been demon- 
strated experimentally to be well separated between elec- 
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tronic recoils and nuclear recoils [nlllj]. This property 
allows these two types of events to be distinguished on 
an event by event basis. The discrimination power is 
proved to be up to 10^ electronic recoils versus one nu- 
clear recoil [ill Ell 20 keV electron-equivalent (ee) 
energy. Such a high pulse shape discrimination (PSD) 
power makes argon a very attractive dark matter target 
for the planned dark matter program at DUSEL. How- 
ever, the existence of ^^Ar and ^^Ar potentially limits 
the detection sensitivity of next generation DUSEL dark 
matter experiments with multi-ton target masses. For 
instance, argon extracted from air results in '^^Ar beta- 
decay at a level of 1 Hz/kg j5j]; a 5-ton detector would 
give a trigger rate of 5000 Hz and a total of 3.2 xlO^'' 
events/year in the energy range of interest (20 keVee to 
60 keVee); this requires PSD to be 3.2 xl0^° in order 
to achieve a sensitivity of 3 x 10"''^ cm^ at 20 keVee 
threshold. Such a high PSD has not been yet demon- 
strated experimentally. Moreover, without depletion of 
■^^Ar, it is unlikely that the data acquisition dead time 
would be manageable for a multi-ton detector and a large 
scale computing facility would be required to perform on- 
line data reduction for a large amount of data from "^^Ar 
decays even for a single-phase detector. In addition, a 
dual-phase detector that drifts electrons from the liquid 
to gaseous would encounter serious difficulties with pile- 
up events because of the milliseconds drift time required. 
With depletion of ^^Ar by a large factor we can not only 
reduce background data, but also make it possible to 
lower the energy threshold to access a larger parame- 
ter space of WIMPs. Therefore, a multi-ton argon-based 
dark matter detector is promising only if natural argon 
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is depleted of '^^Ar. 

A large amount of argon depleted in "^^Ar and '^^Ar 
would be of great benefit to the DUSEL dark matter 
program described above. Because ^^^Ar and ^^Ar are 
produced by cosmic rays in the atmosphere, it is logical 
to consider that underground sources of argon should be 
depleted in '^^Ar and ^^Ar due to the suppression of cos- 
mic rays flux by the overburden. There are two sources 
of underground argon: 1) argon in natural gas and 2) 
argon dissolved in underground water. The former from 
natural gas wells in the US National Helium Reserve was 
found to contain a low level of ^^Ar [l^. The ratio of 
■^^Ar to ^'^Ar (stable isotope) was found to be < 4xl0~^^ 
at a confidence level of 84% [l^ . The latter from under- 
ground water 3200 feet below the surface at Wall, South 
Dakota is the focus of this paper. 

Understanding the production of various argon iso- 
topes underground by stopping negative muon captures 
and neutron interactions is critical to the dating method 
with '^^Ar. In this study, we realize for the first time that 
the ^^Ar production via stopping negative muon capture 
on "^^K is a dominant processes when the depth of un- 
derground resources is smaller than 1800 m.w.e. The 
early work of references showed the '^^Ar produc- 

tion only through (n,p) reaction on '^^K. Neutrons can 
be produced by: 1) cosmic-rays depending on the rock 
composition and depth of the site [15|; and 2) (a, n) neu- 
trons depending on the rock composition and radioac- 
tivity level in the rock. The neutron fluxes and energy 
spectra are essential to the calculation of various isotopes 
of argon and hence the reliability of dating results. In 
this paper, we describe the production mechanisms un- 
derground for argon content in Section |lTl The measured 
rock composition and potassium content in a site as a 
potential argon source is discussed in Section IIIII We 
show the production rate of '^^Ar and "'"Ar in Section ITVl 
Finally, we summarize our conclusions in Section [V] 



1. Stopping negative muon capture on K 

When negatively charged muons are stopped in the 
rock, they either decay through /i^ — > e~j7ejy^ or are 
captured through ijl{A, Z) — )• Z — 1). The capture 

cross section can be estimated using the formula 



where A is atomic mass number of the isotope, Na — 
6.022 X 10^^ mol~^ is Avogadro's constant, v is the speed 
of the stopping negative muons in the medium, r is the 
stopping negative muon capture time in the medium, and 
p is the density of the medium. The lifetime of stopping 
negative muons in potassium is 410 ns [T6| . The capture 
time can be calculated based on 

(2) 

T T„i To 

where r™ (410 ns) is the measured lifetime of muons in 
potassium and tq (2200 ns) is the lifetime of muons in 
vacuum. Using Eqs. (H]) and ([2]), we obtain the cross sec- 
tion of the stopping negative muon capture on potassium 
as a function of the stopping muon energy. 

The production rate of "^^Ar per year per gram of rock 
can be estimated using 

Pi = Na ■ (j) ■ a ■ fneg ' ??, (3) 

where Pi is the production of '^^Ar per year per gram of 
rock, p is the density of potassium in the rock, is the 
stopping muon flux, a is the capture cross section, f„eg 
= 0.44 is the fraction of negative muons [l3|, and 77 is the 
branching ratio of the yield of reaction ^^K{p~ ,1/)^^ Ai . 

2. Neutron-mduced Ar production 



II. THE PRODUCTION MECHANISMS OF 
ARGON ISOTOPES UNDERGROUND 

A. Production of ^^Ar 

Underground production of '^^Ar isotopes occurs pri- 
marily through stopping negative muon captures and 
{n,p) reactions on ■^^K and ^°Ca targets. Both elements 
are in a range of a few percent depending on the chem- 
ical composition of the rock. Neutrons are produced by 
three sources: 1) (a,n) reactions induced by ^38^ g^j^j^ 
^■^^Th decay chains; 2) fission decays induced by ^ssy g^j^j 
232rjij^. gj^^ g-j niuon-induced neutrons. The production 
of '^^Ar through stopping negative muon captures and 
muon-induced neutrons depends strongly on the depth 
of the rock overburden. This results in a depth depen- 
dent production of '^^Ar induced by stopping negative 
muons and muon-induced neutrons. 



In the rock formed thousand years ago, the rates of 
production and decay of '^^Ar should be equal. The pro- 
duction rate per year per gram of rock can be calculated 

by 



P2 



'■ ivr ■ 
^^totC^abs 



(4) 



where (/>„ is the neutron flux in neutrons per year per 
gram of rock, Nj^ and Ntot are the number of potas- 
sium atoms and the total number of atoms in a gram of 
rock, and cr39^(„ and aabs are the exothermic reaction 
■^^K(n,p)'^^Ar cross section and neutron absorption cross 
section. Note that the neutron absorption here means 
that the neutron is absorbed and no further neutron is 
ejected, thus the absorption cross section includes con- 
tributions from all reaction channels which do not eject 
neutrons. 

We calculate reaction cross sections by using the nu- 
clear reaction code TALYS JJj version 0.72. The default 
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input parameters are used in this work. Fig. [T] shows the 
cross section of ^^K(n,p) reaction as a function of neu- 
tron kinetic energy. Also shown in the figure is the cross 
section of *^Ca(n,a) reaction, which can in principle con- 
tribute to '^^Ar production. However, its contribution is 
negligible because of relatively small cross section and 
^^Ca. concentration compared to ^^K. 




Kinetic Energy (l\/1eV) 



FIG. 1: Ar production cross section as a function of neutron 
kinetic energy. 

We plot in Fig. [2] the neutron absorption cross section 
of major atomic components in the rock. The black line 
illustrates the total neutron absorption cross section that 
is the weighted average, which uses the number of atoms 
as the weight factor. 

The number of ^^Ar atoms is given by 



39^,. 



P1+P2 



(5) 



where Pi and P2 are the '^^Ar production rates given by 
Eq.dSl and Eq.(|3]), A = ln2/ii/2 is the decay constant of 
39 Ar. The half-life of ^^Ar is 269 years. 



B. The production rate of ''"Ar 

""^Ar in rock is produced mainly by the decay of ""^K 
through electron capture with a branch ratio of 10.72%. 
The number of ^''Ar atoms produced can be estimated 
according to 



N ^ NoB{l - exp(-Ai)), 



(6) 



where A'o is original atoms at time t = 0, B is the 
branch ratio of electron capture, and A = ln2/ii/2 is the 
decay constant. The half-life h/2 of ^^K is 1.277 x 10^ 
years. Fig. [3] shows the number of ''"Ar atoms produced 
as a function of the rock age. Although the natural abun- 
dance of ''"K is only 0.0117%, the number of ''"Ar pro- 
duced in 65 million years reaches ~ 1.7 x 10'* atoms per 
gram of rock assuming a 2.59% of potassium content in 
the rock. 




Kinetic Energy (IVIeV) 



FIG. 2: The cross section of neutron absorption in different 
nuclei as a function of neutron kinetic energy. The black line 
shows neutron absorption cross section in the rock. 
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FIG. 3: The number of '"'Ar atoms produced in a gram of 
rock as function of the rock age. 



Stopping negative muon capture on *"Ca produces 
*"K, which generates ^''Ar through electron capture de- 
cay. Because the production rate is much smaller than 
the existing in the rock, we ignore this process in our 
calculation of the ratio of '^^Ar to '^^Ar. 



C. Calculation of the ^^Ar/^^Ar ratio as a function 
of depth 

We utilize the fluxes of stopping negative muons and 
muon-induced neutrons as a function of depth derived 
in Ref. [ISJ to calculate the production of argon content 
with equations (|3l4l5l6p . The neutron energy spectrum 
induced by radioactivity obtained in Ref. 19] is used in 
the calculation. 

The common elements found in the Earth's rocks are 
used in the calculation The eight most common 

elements [201 used in the calculation are listed in Ta- 
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ble m The natural radioactivity concentrations of ^ssy 
and ^"^^Th depend on the rock types. We assume the 
common rock in Table |l] is granite. Correspondingly, the 
radioactivity level of 4.7 ppm and 2.0 ppm of ^ssy g^j^j 
^■^^Th are used in the calculation [21j. Using the method 
described in Ref. Ij^, we obtain a yield of 5.5 n/(g-y). 
Subsequently, the '^^Ar production rate is estimated to 
be 7x10"'^ atoms per year per gram of rock. The total 
number of ^^Ar atoms is 2.7 per gram of rock, which is 
independent of depth. 

TABLE I: A list of the eight most common elements found in 
the Earth's rocks. 

Elements Chemical Symbol Mass Weight (%) 



Oxygen O 46.60 

Silicon Si 27.72 

Aluminum Al 8.14 

Iron Fe 5.00 

Calcium Ca 3.63 

Sodium Na 2.83 

Potassium K 2.59 

Magnesium Mg 2.09 



The depth dependence of the "^^Ar production is shown 
in Figure |H Note that the rock chemical composition 
varies from location to location. However, the depen- 
dence of the production of '^^ Ar due to the change in the 
rock chemical composition is much weaker than that of 
the change in depth. 




Depth (km.w.e.) 



FIG. 4: The ^^Ar production as a function of depth. The 
eight common elements found in the Earth's rocks are used 
in the calculation. 

If one assumes that '^^Ar and *°Ar have the same pos- 
sibility to diffuse out of the grains of the rock and to 
be mixed with existing gas in the pore volume, then the 
^^Ar/'^'^Ar ratio in the pore volume can be expressed as 

^ N39Ar X p 

where p is the possibility of '^^Ar and '*°Ar being mixed 
in the pore volume, N^^ is the original number of '^'^Ar 



contained in a certain amount of gas due to the ''"Ar con- 
centration in ambient air, which is related to the pore vol- 
ume, NJ4^ is the number of "^^Ar atoms produced under- 
ground through ^'^K EC decays. Using Eq. ([7]), the ratio 
of the ■'^ Ar to ^"^Ar can be calculated with known natural 
abundance of '^^K in the rock. Given a rock density of 2.7 
g/cm'^, a rock porosity of 20%, and an assumption that 
the pore volume is filled with gas with the ambient air 
concentrations (78%/21%/l% of N2/02/Ar), the ratio of 
■^^Ar to ^°Ar can be determined. The results are shown 
in Figure [5l 




Depth (km.w.e.) 



FIG. 5: The ratio of the ^''Ar to '"'Ar calculated using for- 
mulas and methods described in text. The eight common 
elements found in the Earth's rocks are used in the calcula- 
tion. Shown is a ratio assuming that the produced argon is 
mixed at a 10% level by underground gas. 

Eq.([7]) can be applied to a underground water reser- 
voir, where p is the possibility of '^^Ar and ^°Ar being 
dissolved in water. For a rock density of 2.7 g/cm^, a 
rock porosity of 20%, and an assumption that the pore 
water contains gas with the ambient air concentrations 
(78%/21%/l% of Na/Oa/Ar) at a solubility of 61 mg/L 
for argon, the ratio of '^^ Ar to ^°Ar as a function of depth 
is shown in Figure ID 

As can be seen in Fig. |4l the larger the depth, the 
smaller the production of '^^Ar. On the other hand, at 
the large depth, the production of ^^Ar is dominated by 
stopping negative muon capture and the neutrons from 
natural radioactivity. Since ^''Ca(n, 2p)^^At reaction re- 
quires neutron energy to be greater than 8.3 MeV and the 
majority of the (a, n) neutrons induced by radioactivity 
is less than this reaction threshold, the production of 
^^Ar at large depth underground is dominated by stop- 
ping negative muon capture on '^^K and ■^^K(n,p)'^^Ar. 
Hence, the understanding of "^^K content and radioactiv- 
ity level in the rock surrounding the gas field or water 
reservoir where the argon gas is extracted from is essen- 
tial. 

The underground production of ^'^Ar is primarily from 
^"^K via electron capture decays. Therefore, an accurate 
knowledge of the potassium content is critical. The un- 
derground production of '^^Ar is negligible. 
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FIG. 6: The ratio of the ^^Ar to ""Ar calculated using for- 
mulas and methods described in text. The eight common 
elements found in the Earth's rocks are used in the calcula- 
tion. Shown is a ratio assuming that the produced argon is 
dissolved at a 10% level by underground water. 



III. EXPERIMENTAL RESULTS OF 
POTASSIUM CONTENT 

A. Underground water Site 

The community of Wall, South Dakota pumps water 
using seven vifells that extract water from the Fall River 
aquifer in an old rock formation 3200 ft below the sur- 
face. The average water consumption is 200 gallons per 
minute. The water temperature underground is about 
42°C. The expected argon gas extraction is about 2.3 kg 
per day assuming 70% efficiency with 61 mg/L solubil- 
ity. However, the level of '^^Ar in argon gas extracted 
from well water is contingent upon the understanding of 
the rock composition especially the potassium content in 
the rock. Two different samples obtained from the South 
Dakota Geological Survey were evaluated. The first sam- 
ple was a cutting sample from the 3200 ft level in the well 
that was provided by the well drilling company. The sec- 
ond sample was a core sample taken from the Fall River 
Aquifer at a distance of about 40 miles from Wall. 



B. Method 

The samples were broken down using a mortar and 
pestle and then passed through a 100 mesh sieve. The 
samples were then dried in an oven set at 100°C to re- 
move the water. After drying the samples were consis- 
tently stored in a desiccator to prevent re-hydration. Ini- 
tially the amount of free K (K not in the soil matrix) was 
found using the following solution preparation and uti- 
lizing the atomic absorption spectrometer. The solution 
preparation consisted of 2 g of the sample, 40 ml of water 
(Note that all water used was nanopure water that had 
been purified to 18 Mfl cm with a SYBRON/Barnsted 
Nanopure System), 15.5 ml of concentrated hydrochloric 



acid (HCl), and 2 ml of concentrated nitric acid (HNO3 
). One important observation was that the solution did 
not boil when the acids were added indicating a lack 
of carbonates. The solution was heated with stirring, 
and finally filtered to remove the remaining soil particles. 
The filtered solution was placed in a 250 ml volumetric 
fiask and diluted with water to the 250 ml mark. An 
atomic absorption analysis was done on this sample us- 
ing a Thermo Jarrell Ash atomic absorption spectrometer 
which contained a Smith-Hieftje 11 spectrophotometer. 
The light source was a hollow cathode lamp produced 
by the Thermo Jarrell Ash Corporation which emitted 
at 766.5 nm. The bandpass was set at 1.0 nm and the 
lamp current was set to 5 mA. Two runs were done with 
the instrument for each sample. The measurement of the 
total K in the soil samples was measured under the as- 
sumption that most of the K was tied up in a Si02 com- 
pound. This was concluded from the X-ray diffraction 
analysis conducted on the sample. The sample prepara- 
tion was obtained from Buck Scientific (2^ . One gram 
of the soil sample was mixed with 5 grams of a 2:1 fiux 
of lithium carbonate (LiCOa) and zinc oxide (ZnO). The 
sample and flux were placed in a platinum crucible and 
fused in a muffie furnace that was heated to 950° C for 
20 minutes. The melt was then allowed to cool and was 
re-dissolved using 20 ml of HCl and 20 ml of water. This 
solution was filtered to remove the soil particles that had 
not been re-dissolved. A series of standards that ranged 
from 0.01 ppm to 5.0 ppm were prepared and analyzed 
with the atomic absorption spectrophotometer yielding 
a calibration curve of y = (0.0562 ± 0.0009)a; - (0.0034± 
0.0021) with a correlation coefficient of 0.9994. Potas- 
sium samples were diluted to fall within the calibration 
curve. Once again two runs were done and then averaged. 
Three different samples from the core sample and from 
the cutting sample were treated in this manner. Finally, 
we obtain a potassium content of 18500 ± 300 ppm. 



IV. PRODUCTION RATE OF ARGON 
ISOTOPES 

In this particular site in South Dakota at 3200- ft below 
the surface, the stopping muon flux is calculated to be 
3.4xl0"^cm"2s"^ Utilizing Eq.®, we obtain a ^'^Ar 
production of 0.0026 per year per gram of rock. 

The natural radioactivity of ^"^^Th and ^ssy jg 
lyzed to be around 0.6 ppm and 1.9 ppm, separately. The 
rock is sandstone with a chemical composition described 
in Table HH 

The (a, n) neutron yield in the rock is calculated [l^ 
to be 0.25 n/(g-ppm-y) for ^^^Th and 0.65 n/(g-ppm-y) 
for ^^^U. The fission decay of ^asy induced neutron yield 
is 0.52 n/(g-ppm-y). The total neutron yield induced 
by radioactivity is 2.4 n/(g-y). Using Eq.(|3]), the '^^Ar 
production rate is estimated to be ^3x10""^ atoms per 
year per gram of rock. 

The '^^Ar produced by the muon-induced neutrons at 
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TABLE II: The chemical composition of the sandstone 3200 
ft below the surface at Wall, South Dakota. 




this depth is negligible. The total number of '^^Ar atoms 
given by Eq.([S]) is estimated to be 2.17 atoms per gram 
of rock. 

If we assume that '^^Ar and '*"Ar have the same pos- 
sibility to diffuse out of the grains of the rock and to be 
dissolved in the pore water, then the ■^^Ar/^'^Ar ratio in 
water can be obtained by applying Eq. ([7]), 



NAr + 1.2 X 1014p ^ ^ 

where p is the possibility of '^^Ar and '"^Ar being dis- 
solved in water, and NAr is the original number of '^'^Ar 
contained in a certain amount of water, which is related 
to the pore volume. If we assume that the porosity is 
20% (this is typical for sandstone ) , then the water vol- 
ume corresponding to a gram of rock with the measured 
density of 2.2 g/cm^ is 9.1 x 10~^ cm''. Given a total 
saturation solubility of 125 mg/L for N2/02/Ar at 20°C 
assuming the gas has ambient air concentrations (nitro- 
gen/oxygen/argon of 78%/21%/l%) with different solu- 
bility of N2 (20 mg/L)/02 (44 mg/L)/ Ar (61 mg/L), the 
dissolved argon in water is 61 mg/L x 0.01 / (20 mg/L 
X 0.78 + 44 mg/L x 0.21 -I- 61 mg/L x 0.01). This re- 
sults in a 2.4 % of argon content ( 3 mg/L) in the total 
gas content in the water. Therefore, NAr — 4.1 x 10^^. 
The term of 1.2x10^^ is the number of ''"Ar atoms from 
''"K decay with a total concentration of 18500 ppm in a 
~65 million- year rock formation. Fig. [7] shows the ra- 
tio of '^^Ar to '*°Ar atoms contained in water as a func- 
tion of p. Thus, a depletion factor of 1.6 to 155 can be 
achieved for the argon extracted from this site when p 
varies from 0.01 to 1.0. Considering saturation solubility 
decreases when temperature increases, the dissolved pos- 
sibility of '^^Ar and ^"Ar can be as small as 1% for the 
argon content produced in the underground environment 
where the temperature is higher. This would make un- 
derground water reservoir a valuable source for depleted 
argon that can be used for ton-scale dark matter exper- 
iments. If one assumes that 30% [3 of Ar produced will 
be dissolved in water, the ratio between ^^Ar and ^'^Ar 
in water is ~ 1.6 x 10~^^, which is a factor of ^-^5 lower 
than that in atmosphere. 



FIG. 7: The ratio of Ar to Ar atoms contained in water 
as a function of the dissolved probability (p). 



V. CONCLUSION 



We have evaluated that the '^^Ar production under- 
ground through (n,p) reactions induced by neutrons and 
stopping negative muon capture. We demonstrate that 
argon depleted in '^^ can be obtained only if the depth of 
the underground resources is greater than 500 m.w.e. be- 
low the surface. Stopping negative muon capture on "'^K 
dominates over radiogenic production at depths of less 
than 2000 m.w.e., and that production by muon-induced 
neutrons is subdominant at any depth. Depending on the 
rock porosity and argon diffusion coefficient, the upper 
limit of the concentration of ^^Ar in the underground 
water at Wall of South Dakota is estimated to be in a 
range of a factor of 1.6 to 155 less than the ^^Ar concen- 
tration in the atmosphere. Naturally occurring EC 
decay dominates the content of ^''Ar underground. The 
ratio of ■'^Ar to "'"Ar depends largely on the following: 
1) the potassium content in the rock; 2) the natural ra- 
dioactivity of ^^^Th and ^ssy ^j^g j-Qck; 3) the rock 
porosity; 4) argon gas diffusion coefficient underground; 
and 5) argon gas solubility in the water underground. 
We estimate that the production of ^^Ar is negligible at 
this underground site. The calculation tools presented 
in this paper are also critical to the dating method with 
39 Ar IIQ m, [23. We conclude that the ratio of ^^Ar 
to ''"Ar varies from location to location depending on 
the content and the radioactivity level in the rock at 
large depth. Extraction of argon gas from underground 
sources must: 1) evaluate the depth of the site, the age 
of the rock, the rock composition, and radioactivity; 2) 
study the porosity of the rock and argon gas diffusion co- 
efScient; and 3) understand the solubility of argon gas as 
a function of pressure and temperature. The algorithm 
described in this paper is critical to understanding the 
discrepancy in the dating methods with the ratio of ''"Ar 
to 4"Ar [11 [M^. 
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